London penetration depth and strong pair-breaking in iron-based superconductors 
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The low temperature variation of the London penetration depth for a number of iron-pnictide and 
iron-chalcogenide superconductors is nearly quadratic, AA(T) = j3T n with n ~ 2. The coefficient 
in this dependence shows a robust scaling, (3 oc 1/T% across different families of these materials. 
We associate the scaling with a strong pair-breaking. The same mechanism have recently been 
suggested to explain the scalings of the specific heat jump, AC oc TJ? and of the slopes of 
the upper critical field, dHcijdT oc T c in these materials Q]. This suggests that thermodynamic 
and electromagnetic properties of the iron-based superconductors can be described within a strong 
pair-breaking scenario. 

PACS numbers: 74.20.-z, 74.20.Rp 



Due to the unique electronic structure and, most likely, 
unconventional pairing mechanism, iron-based supercon- 
ductors exhibit a number of uncommon properties. It has 
recently been reported [l[ that across the whole family of 
iron-pnictides the specific heat jump, AC, at the critical 
temperature T c shows an extraordinary scaling AC oc Tj?, 
whereas in conventional s-wave materials AC oc T c . Ac- 
cording to Ref.0, this unusual scaling is caused by a 
strong pair-breaking in materials with anisotropic order 
parameters; both transport and magnetic scattering in 
such materials suppress T c and there are plenty of rea- 
sons for magnetic pair breaking in iron-based supercon- 
ductors. Another consequence of this model, proportion- 
ality of slopes of the upper critical field [dH c2 /dT]T a to 
T c , has also been shown to hold for the data available Q. 
In this work we show that the same idea can be applied 
to the low temperature behavior of the London penetra- 
tion, A A = A(T) — A(0), where the pair-breaking results 



AA oc T 2 /T 3 



(1) 



Despite some initial disagreements in experimental re- 
ports, most precision measurements of the in-plane Lon- 
don penetration depth of iron-based superconductors had 
found the power-law behavior, AA(T) oc T n with n ~ 2 
for some compounds n 
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1 is claimed 

Commonly, a non-exponential behavior is taken as ev- 
idence of unconventional order parameter, possibly hav- 
ing a nodal gap structure [1, 0, 0, EH • However, such a 
direct correspondence between the nodes and the expo- 
nent n should exist only in clean materials. As a rule, 
scattering breaks this elegant connection. E.g., for d- 
wave superconductors, the linear low T dependence of A 
in the clean case changes to T 2 in the presence of mod- 
erate scattering 14j. In fact, connection between the 



power-law behavior of AA(T) and scattering in pnictides 
had been suggested @, 0,1, El- The symmetry of the 
order parameter A in multi-band iron-pnictides is not yet 
determined with certainty, however, many favor the s± 
structure 15, T^. The Fermi surface (FS) average of the 



order parameter in this model (A) <C A max . We then 
expect the penetration depth to behave like a "dirty" d- 
wave, i.e., to show the low-temperature variation oc T 2 . 

The samples were plate-like single crystals with typical 
dimensions 1 x 1 x 0.2 mm 3 . Details of sam ple syn thesis 
and characterization can be found elsewhere 17H19I] . The 
penetration depth measurements were performed with a 
sclf-rcsonating tunnel diode oscillator. Diamagnetic re- 
sponse of the sample causes shift of the resonant fre- 
quency, A/ = —Gx(T), where x(T) is magnetic sus- 
ceptibility determined by A(T) in the Meissner state, 
-4ttx(T) = [1- (A/-R)tanh(-R/A)]. The calibration con- 
stant G = foV s /2V c (l — N) is measured directly by ex- 
tracting the sample from the coil at the lowest temper- 
ature. Here /o ~ 14 MHz is the empty resonator fre- 
quency, V s and V c are the sample and coil volumes, and N 
is the demagnetization factor. Details of measurements 
and of data analysis are described elsewhere 2(| 21 1 

Figure [3a) shows the linear behavior of AA versus 
(T/T c ) 2 for T < T c /3 in iron-based compounds with T c 
varying from « 12 to 23 K; the data are from Refs. 0- 
7, and 22]. The exponent n in AA oc T™ extracted by 
fitting the low temperature data is shown for six com- 
pounds in Fig.QJb). We see that AA(T) oc T 2 holds for 
the (AE)(Fei_ :c TM :c ) 2 As 2 ("122"), (RE^eAs^i^F,) 
("1111"), and FeTei-^Se^ ("11") families; here AE 
stands for an alkali earth element, TM for a transition 
metal, RE for a rare earth. Thus, the four lines shown in 
Fig. [Ha) are not merely for different doping levels of the 
same compound, but rather they belong to four different 
families of the iron-based materials. This universal be- 
havior has prompted us to look for a universal cause; we 
offer below a strong pair-breaking as such a cause. 

The theoretical tool we employ, the quasiclassical ver- 
sion of the weak-coupling Gor'kov theory, holds for a gen- 
eral anisotropic Fermi (F) surface and for any gap sym- 
metry [23I ] . The formalism in the form convenient for our 
purpose is outlined in Ref. we refer readers to this work 
for details. The theory is formulated in terms of functions 
f(r,kp,oj), f + , and g which originate from Gor'kov's 
Green's functions and are normalized by g 2 + ff + = 1; 
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FIG. 1. (Color online) (a) AA versus (T/T c ) 2 

for Ba(Feo.942Coo.o58)2As2 marked by (1), 
Ba(Feo.g4iNi .o59)2As 2 (2), Fe1.001Seo.367Teo.632 (3), and 
LaFcAsOo.gFo.i (4). Inset: (a) AA in the full temperature 
range, (b) Fitted exponent n in A A oc T n . 



the Matsubara frequencies are to = itT(2v + 1) with an 
integer v and % = kg = 1- The order parameter is taken 
in the form A(r,k F ) = ^(r,T)Q(k F ) where Q(k F ) de- 
scribes the variation of A along the F-surface and is con- 
veniently normalized so that the average over the whole 
F-surface (ft 2 ) = 1. Hence, the model is a BCS-type 
weak-coupling approach providing a qualitative descrip- 
tion at best. 

The scattering in the Born approximation is character- 
ized by two scattering times, the transport r responsible 
for the normal conductivity and r m for processes break- 
ing the time reversal symmetry (e.g., spin- flip): 

l/r± = 1/r ± l/r m . (2) 

Commonly, two dimensionless parameters are used: 

p = 1/2ttT c t and p m = l/2nT c T m , (3) 

or equivalently p± = p ± p m . This is of course a gross 
simplification. For multi-band F-surfaces one may need 
more parameters for various intra- and inter-band pro- 
cesses, which are hardly controllable and their number 
is too large for a useful theory. Our model is amenable 
for analytic work and may prove helpful, the simplifying 
assumptions notwithstanding. 



It is well-known that the formal scheme of the semi- 
nal Abrikosov-Gor'kov (AG) work on magnetic impuri- 
ties Q applies to various situations with different pair- 
breaking causes, not necessarily the AG spin-flip scat- 
tering 25]. In each particular situation, the parameter 
p m must be properly defined. Here, without specifying 
the pair breaking mechanism, we apply the AG approach 
to show that the pair-breaking accounts for our data on 
the low temperature X(T) along with the earlier reported 
behavior of H C 2 slopes at T c and of the quite unusual de- 
pendence of the specific heat jump on T c . 

Evaluation of A(T;r, r m ) for arbitrary r's and arbi- 
trary anisotropy of A is difficult analytically. However, 
for a strong T c suppression, the problem is manageable. 
Within the microscopic theory, penetration of weak mag- 
netic fields into superconductors is evaluated by first solv- 
ing for the unperturbed zero-field state and then treat- 
ing small fields as perturbations. It was shown by AG 
24| that for strong pair-breaking the formalism for the 
derivation of the Ginzburg-Landau equations near T c ap- 
plies at all temperatures. Within the Eilenberger ap- 
proach this means that /« 1 and - // + /2 at all 
temperatures. The calculation then proceeds in a manner 
similar to that near T c . 

Within a two-band model for iron-based materials, the 
order parameter is believed to have a ±s structure [l5j ]. 
so that (A) <C |A maa; | The problem is considerably 
simplifies if one assumes (A) = 0; we use this assumption 
and expect the model to hold at least qualitatively. In 
the zero-field state, we look for solutions of Eilenberger 
equations as / = + / (2) + ... where /W ~ A, /( 2 ) ~ 
A 2 , etc. The Eilenberger equation for / then yields [2|]: 



fo 



A 



A / (A 2 ) 
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1/2T+. 
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0(A 5 ), (4) 



One can see that even at low 
t+T c ~ l/p+ <C 1 because for 



where cj+ = lu - 
temperatures fo 
strong pair-breaking T c — > 0. This is a quasiclassical 
justification for the AG statement that / <C 1 at all T's. 

The T dependence of A (or 'J) is obtained with the 
help of the self-consistency equation (or the "gap equa- 
tion"). For a strong pair- breaking, this equation takes 
the form 0]: 



E 
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flf 



(5) 



Substituting here / of Eq. ^ , we obtain the order pa- 
rameter in the field-free state: 



27r 2 (r 2 -r 2 ) 

3(f2 4 ) - 2 



(6) 



this reduces to the AG form for Q = 1. 

We can now consider the response to a small current 



j = -47r|e|AT(0)T Im ^ (vg 



w>0 



(7) 
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N(0) is the density of states at the F- level per one spin. 
Weak supercurrents leave the order parameter modulus 
unchanged, but cause the condensate to acquire an over- 
all phase 9{r). We then look for perturbed solutions as: 



/ + = (/o + tf)e* 



(fo + fi)e w , 
9 = 9a + 9i, 



(8) 



where the subscript 1 denotes small corrections to the 
uniform state /o,<?o- I n the London limit, the only coor- 
dinate dependence is that of the phase 9, i.e., fi,gi are 
r independent. The Eilenberger equations provide the 
corrections among which we need only g\ : 



9i 



2lu + 



iA 2 vP 

2cui 



(9) 



see 0. Here P = V6 + 2vrA/0 o = 2na/<f>o with the 
"gauge invariant vector potential" a. 

We now substitute go + g\ in Eq. and compare the 
result with Airji/c — — (A 2 )7L Ofc to obtain: 



87T 2 e 2 A^(0)T c 



ViVkO,' 



(10) 



OJ>0 



The sum here is expressed in terms of the polygamma 
function: 



uj>0 



1_ 1 I 

xd ~ 167T 3 T 3 ^ I 2t + 2 



^ <"> 



where p+ ^> 1 has been used. 
Eq. ||(5J), one obtains: 



Taking into account 



(A 2 



167r 3 e 2 7V(0)fc 
c 2 ft 2 (3(0 4 ) 



2 T 2 

' ^^ fe Q 2 )(T c 2 -T 2 ) (12) 



in common units. It is now easy to obtain the low T be- 
havior of AA b = X a b(T) — A a fc(0) for a uniaxial material: 



T 2 



eh 



rj = 



3(f> 4 } 



8nk B T + I/ 7re 2 iV(0)(w 2 fi 2 ) 



.(13) 



We stress that r + here is close to the critical value for 
which T c — > 0. One readily obtains for T = 0, 



A ab (0) = 2 V /T C 



(14) 



Note: Eqs. fFSfr and (H]) are derived for (0) « 0. One 
can show that they hold for (f2) ^ as well with, how- 
ever, different coefficient rj. We do not provide here this 
cumbersome calculation. 

To examine the predicted scaling behavior, the factor 
(3 in AA = j3T 2 was obtained by fitting the low temper- 
ature AA for the same 122, 1111 and 11 compounds of 
Fig. Q] with (3 being the only fitting parameter. The /3's 
are plotted in the main panel of Fig. [5] versus T c . The 
error bars on this graph reflect the fact that each sam- 
ple studied has a certain transition width. The inset of 
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FIG. 2. (Color online) The factor /3 obtained in fitting of data 
to AA = /3T 2 plotted versus T c on a log-log scale. The solid 
line is a fit to (3 = rj/T^, motivated by Eq. (|13[l for a strong 
pair-breaking. 



Fig. [2] shows the convention adopted here for T c determi- 
nation. The uncertainty of T c is the dominant source of 
error in determination of (3. According to the strong pair 
breaking scenario, (3 — rj/T^. To compare experiment 
with theory, (3 is plotted on a log-log scale in the main 
frame of Fig.[2] along with the line (3 = (8.8 ± 1.0)/T 3 
obtained by fitting the data. Moreover, by substituting 
v ~ 10 7 cm/s and N(0) ~ 10 33 erg^cm -3 in Eq. (O 
we roughly estimate r + ~3x 10~ 14 s; this value corre- 
sponds to parameter p + w 5 for T c = 40 K and to larger 
values for lower T c 's, an observation consistent with the 
major model assumption of p + 3> 1. The degree to which 
the experimental values follow the theory is remarkable, 
a substantial scatter of the data points notwithstanding. 

The scalings of Eqs. (fT5j) and (TTJ} are obtained for a 
strong pair-breaking materials with the order parame- 
ter obeying (f2) w 0. These conditions are likely to 
be satisfied in underdoped high-T c cuprates since under- 
doped materials are clearly disordered and the d-wave 
order parameter is suppressed by any scattering. In- 



deed, the surface resistance [26| and optical data [27J for 
YBa2Cu306+£c samples with T c varying from 3 to 17 K 
show 1/A 2 fc (0) oc T 2 in agreement with Eq. (JUJ). This 
behavior differs from "Uemura scaling" 1/A 2 b (0) oc T c 

To our knowledge there is not yet sufficient data on 
A(0) for the iron-based materials to verify the scaling 
(fl4)) . Similarly, we are not aware of a data set to check 
the scaling H c \ tx T 2 which follows from Eq. (fT4|) . 

We would like to stress that the penetration depth 
scalings discussed here as well as those for the the spe- 
cific heat jump and for the slopes of H c2 (T) described 
in Ref. [2| are approximate by design since their deriva- 
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tion involves a number of simplifying assumptions. Still 
they are robust in showing that the pair-breaking is an 
important factor in superconductivity of iron-pnictidcs. 

Many questions still remain; for example, why the Co 
doped 122 compounds deviate substantially from the gen- 
eral scaling behavior shown in Fig.[5J see also Ref . Q • An- 
other problem to address is how to reconcile the strong 
pair-breaking, which in the isotropic case leads to gapless 
superconductivity [24|, with the in-plane thermal con- 
ductivity data showing k(0) = 2tJ 3(|. At this point, 
we can say that (a) the strong pair-breaking model for 
anisotropic order parameters states that the total density 
of states N(e) integrated over all pockets of the Fermi 
surface is finite at zero energy [2|; this does not exclude 
a possibility that N = for some parts on the Fermi sur- 
face. And (b): in this work we are interested in the su- 
perfhhd density oc 1/A 2 which depends only on the Fermi 
surface average (A) so that our results are less sensitive 
to the A behavior on a particular set of directions (e.g., 
those in the ab plane). The same qualitative argument 
shows that our scalings do not contradict the in-plane 
ARPES data [H|. 

To conclude, analysis of the low-temperature behavior 
of the London penetration depth shows that a strong 
pair-breaking is likely to be responsible for the nearly 
universal temperature dependence AA Q b oc T 2 /Tj?, along 
with earlier reported AC oc T~ 3 and [dH C 2/ 'o]T]t c oc T c , 
in nearly all iron-based superconductors. 
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